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Abstract 

Nucleon properties are investigated in background electric fields. As the magnetic moments 
of baryons affect their relativistic propagation in constant electric fields, electric polarizabilities 
cannot be determined without knowledge of magnetic moments. This is analogous to the experi- 
mental situation, for which determination of polarizabilities from the Compton amplitude requires 
subtraction of Born terms. With the background field method, we devise combinations of nucleon 
correlation functions in constant electric fields that isolate magnetic moments and electric polariz- 
abilities. Using an ensemble of anisotropic gauge configurations with dynamical clover fermions, we 
demonstrate how both observables can be determined from lattice QCD simulations in background 
electric fields. We obtain results for the neutron and proton, however, our study is currently lim- 
ited to electrically neutral sea quarks. The value we extract for the nucleon isovector magnetic 
moment is comparable to those obtained from measuring lattice three-point functions at similar 
pion masses. 
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I. INTRODUCTION 



Understanding low-energy properties of hadrons directly from QCD remains a challeng- 
ing endeavor. In this low-energy regime, quark and gluon interactions must be treated non- 
perturbatively, ultimately resulting in their confinement into hadrons. After three decades 
of dedicated work, lattice QCD has evolved into a tool to address quantitatively the non- 
perturbative dynamics underlying hadrons and their interactions, see |l| for an overview. 
Electromagnetic moments and multipole polarizabilities are low-energy properties of hadrons 
with transparent physical meaning. These properties characterize the distribution of charge 
and magnetism within a hadron, and the response of the charge and magnetism distributions 
to external fields, respectively. Low-energy properties of hadrons can be described using an 
effective theory of QCD, based upon treating pseudoscalar mesons as the Goldstone modes 
arising from spontaneous chiral symmetry breaking. A picture of hadrons emerges from chi- 
ral dynamics: that of a hadronic core surrounded by a pseudoscalar meson cloud. In part, 
the electromagnetic properties of hadrons encode the distribution of charged mesons, and 
the stiffness of the charged meson cloud. Chiral dynamics consequently makes predictions 
for the form of electromagnetic observables. Confirming these predictions both experimen- 
tally and from the lattice will be a milestone in our understanding of non-perturbative QCD 
dynamics. 

Computation of hadronic electromagnetic properties using lattice QCD can be accom- 
plished in at least two different ways. The current insertion method, see e.g. can be 
used to determine hadronic matrix elements of the electromagnetic current. This method is 
ideal for the computation of electromagnetic form factors, but is limited in the extraction 
of multipole moments due to the available lattice momentum, which, for periodic boundary 
conditions, is quantized in units of 27r/L, where L is the size of the lattice.^ For current lat- 
tice sizes, the determination of multipole moments relies on a long extrapolation to vanishing 
momentum transfer. For multipole polarizabilities, the temporal extent of current lattices 
also makes direct computation of the Compton scattering tensor infeasible. Were lattices 
long enough to allow the computation of matrix elements with two current insertions, the 
extraction of polarizabilities would still require a long extrapolation to zero momentum. Al- 
ternately, hadronic properties can be determined using the background field method [sl-H]- 
With this method, one determines lattice two-point functions in the presence of classical 
external fields. Observables are then determined from the variation of these hadronic corre- 
lators with the strength of the external field. 

In the current context, we focus our attention on properties of spin-half baryons in 
external electric fields. Neutral hadrons in electric fields have been investigated with lattice 
QCD using the quenched approximation at pion masses greater than 500 MeV js], lof. A recent 
calculation has explored neutron properties at lower pion masses with electrically neutral 
sea quarks [lo|. There has also been a fully dynamical calculation for the neutron in an 
electric field using a pion mass of 760 MeV . In this work, we treat the baryon spin in 
a relativistic manner, and thereby demonstrate how to determine both nucleon magnetic 
moments and electric polarizabilities using lattice QCD in background electric fields. In 
essence, our method allows one to perform the subtraction of the nucleon pole term aris- 



For isovector form factors, the restriction to quantized momentum transfer can be lifted by imposing 
isospin twisted boundary conditions on the quark fields 0] ■ The only known method to handle the 
isoscalar contribution is to increase the lattice volume. 
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ing from anomalous magnetic couplings. The resulting polarizabilities consequently have 
the correct physical interpretation, and correspond to those extracted from experiment, as 
well as those derived from chiral perturbation theory, for example. We consider both the 
neutron and proton in this study. For the latter, we use the relativistic generalization of 
the method proposed in which relies on matching QCD correlators onto those derived 
from single-hadron effective actions. We have also employed this method recently for pseu- 
doscalar meson electric polarizabilities [l3|. A salient feature of these computations is that 
they utilize a periodic lattice action with everywhere constant electric fields, and thereby 
eliminate difficulties arising from Dirichlet boundary conditions used in previous studies.^ 
Our calculations of nucleon magnetic moments and electric polarizabilities include effects 
from dynamical quarks, however, they are restricted to electrically neutral sea quarks. 

We organize our presentation in the following manner. First in Sect. [Tll we analytically 
determine the form of baryon correlation functions in external electric fields. We specialize 
to the case of a uniform electric field, and derive results for both neutral and charged 
baryons. A key observation of this section is that baryon electric polarizabilities cannot 
be determined without knowledge of their magnetic moments. This is analogous to the 
experimental situation, where Born terms must be subtracted to extract polarizabilities 
from Compton scattering. Appendix |A] is concerned with the physics underlying the Born 
subtraction. In Sect. IIIIl we provide the pertinent details of our lattice computations, 
and implementation of the background field. In Sect. llVt we present our analysis of nucleon 
correlation functions calculated in background electric fields using lattice QCD. For both the 
neutron and proton, we demonstrate that the measured correlation functions agree in form 
with the analytic expectations from the hadronic theory, and that magnetic moments and 
electric polarizabilities can be extracted from data. Appendix |B] is devoted to the analysis of 
unpolarized neutron correlators, from which consistent results are obtained. (These results, 
however, do not permit the determination of the electric polarizability — only a combination 
of the polarizability and the square of the magnetic moment. The latter contribution arises 
from Born-level couplings.) A brief conclusion in Sect. |V]ends our work. 



II. SPIN-HALF CORRELATION FUNCTIONS 



To extract properties of nucleons in background electric fields, we must first understand 
the expected behavior of their two-point correlation functions. In this section, we determine 
baryon two-point functions using the single-hadron effective action that arises from QCD in 
the ultra- low energy limit. The functional forms deduced for these two-point functions can 
then be utilized to fit baryon correlators computed with lattice QCD. From these fits, one 
can deduce hadronic parameters, such as the magnetic moments and electric polarizabilities. 

To arrive at a uniform electric field of the form £ = Sz, we use the Euclidean space 
vector potential 

= (0,0,-^X4,0). (1) 
While there are other gauge equivalent choices, we find Eq. ([1]) particularly useful.^ The 



Periodic actions with everywhere constant magnetic fields have also been employed recently to study 



15, 



16|. 



magnetic moments of hadrons IJj, and modification of the QCD vacuum 
^ On a torus, many of the gauge equivalent choices in infinite volume are no longer equivalent, but differ 
by their holonomy. External fields with non- vanishing holonomy lead to new interactions that are finite 
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analytic continuation, 8 — )■ —zEm, is needed to recover Minkowski space results. As our 
interests lie only with quantities perturbative in the strength of the field, this analytic con- 
tinuation can be performed trivially, see jl^l- A Euclidean formulation is natural from the 
point of view of lattice gauge theory simulations; moreover, the Euclidean formulation re- 



moves instabilities due to non-perturbative effects, i.e. the Schwinger mechanism j2l|. With 



the vector potential specified, we can determine the baryon two-point functions. As neutral 
and charged baryons propagate differently in electric fields, we handle each separately. 



A. Neutral Spin-Half Baryons 

We consider first the case of a neutral spin-half particle of mass M described by the 
field ip{x). The Euclidean space correlation function in the hadronic theory we denote by 
Gai3{x4^,S), which is given by 



G^p{x4,£) = J dx{0\Mx)^^{Om£, (2) 

where the subscript denotes that the correlation function is calculated in the background elec- 
tric field. Integrating over all space projects the correlator onto vanishing three-momentum, 
which is a good quantum number. For lattice QCD with spatially periodic boundary condi- 
tions, a sum over lattice sites accomplishes the same thing. 

The energy of the neutral particle, E{S), depends on the strength of the electric field. 
For weak fields, the energy has the expansion 

E{S) = M+U7iaES^ + ..., (3) 

where a^; is the electric polarizability, and the ellipsis denotes higher-order terms in even 
powers of the field. The quadratic Stark shift is positive due to our Euclidean space treat- 
ment. The magnetic moment, fi, is also important, and this coupling is entirely anomalous. 
The single-hadron effective action for a neutral spin-half particle takes the form 



Se = d'^x ijj{x) 



^(x). (4) 



This action describes the dynamics of the neutron in the ultra-low energy limit of QCD. The 
electromagnetic field strength tensor is F^i, = d^A^, — dyA^j,. For a background electric field, 
(^iiuF^y = 2K ■ S, where K = i^'j^ is the generator of boosts in the spin-half representation 
of the Lorentz group. The magnetic moment coupling has been written as £^-dependent. In 
small fields, fJ^iS) has a perturbative expansion in even powers of the field, and satisfies the 
zero-field limit, fi{0) = /i. Using the effective action in Eq. (jlj) to determine the unpolarized 
two-point function, we arrive at 

Tr[G'(x4,^)] = Z{S)exp[-x,E,4S)], (5) 



As the current study is restricted to one lattice volume, we postpone the 



volume artifacts 
investigation of finite volume effects to future work. 
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where the effective energy, Ees{£), depends on the magnetic moment, and is given by 



= M + i£=(4™,-J^)+.... (6) 

In the second hne, we have retained terms in the effective energy only up to second order in 
the electric field. There are two distinct contributions at that order. The first is the expected 
shift in energy due to the polarizability, while the second is the analogue of the nucleon pole 
term arising from Born-level couplings to the magnetic moment. Appendix |A] details the 
physics underlying this Born-like contribution. While the magnetic moment contribution 
arises from a relativistic effect, the resulting shift in baryon energy occurs at the same order 
as that due to the polarizability.^ Appendix [B] is devoted to the extraction of Ees{S), using 
unpolarized neutron correlation functions. The correlator in Eq. ([5]), however, does not 
allow access to the electric polarizability without knowledge of the magnetic moment. A 
background field analogue of the Born subtraction is needed. 

To extract both the magnetic moment and electric polarizability, we use the boost 
projection operators 

V± = ^{1±K,), (7) 

where is the boost operator is the z-direction. The boost-projected correlation functions 
are given by 

G±{x^,S) = Tt[V±G{x,,£)] (8) 

With the additional electric field dependence present in the amplitude, one can separate the 
electric polarizability from the magnetic moment by simultaneously analyzing both boost- 
projected correlators. The method we employ to accomplish this will be detailed below. 



B. Charged Spin-Half Baryons 



Consider now a spin-half baryon with charge Q. The magnetic moment, n, is a sum of 
two terms, fi = Q + jl. The piece proportional to the charge is the Dirac magnetic moment, 
while that denoted by jl is the anomalous magnetic moment. Including terms relevant for 
a uniform external field, the relativistic single-particle action for a charged baryon has the 
form 

Se = I d*x^{x) p+ E{£) - ^^^^ 



ip{x) 



(10) 



^ In the chiral limit, the contribution from the magnetic moment term is suppressed relative to the electric 
polarizability by a factor of m^r/Af. This suppression, however, owes to the singular behavior of the 
polarizability in that limit, namely ue ^ I/tjIt^. Because our lattice pion mass is larger than physical, we 
will make no assumption about the dominance of the polarizability near the chiral limit. 
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where the electromagnetically gauge covariant derivative is = dfj, + iQA^j^, and higher- 
order terms in the field strength appear parametrically in the £^-dependent couphngs in 
Eq. (fTOj) . The parameter E{S) is the charged particle's rest energy, which has the weak field 
expansion in Eq. ([3]). The anomalous magnetic moment coupling, p^{S), has a weak field 
expansion in even powers of the field, and satisfies the zero-field relation, fl{0) = jl. The 
action in Eq. f llOp describes the proton in the ultra-low energy regime of QCD. 

As with the neutral baryons, it is beneficial to consider boost-projected correlation 
functions. For charged baryons, these have the form 



G±{x4^,£) 



Z{£) [l ± D (x4, E,^{£f T Q£, £) 



where the function D{x4,E'^,£) is the relativistic propagator function 20 



(11) 



poo 

D{xi,E'^,£) = J ds^ 



Q£ 



27r smh{Q£s] 



■ exp 



~Q£xlcoth{Q£s) - ^E'^s 



:i2) 



When the charge is set to zero, we recover the neutral baryon correlation functions in Eq. ((91). 
As one can see, to determine charged baryon electric polarizabilities, we must also deduce 
their anomalous magnetic moments. This can be achieved utilizing both boost-projected 
correlation functions. These fits are more complicated than for neutral hadrons; however, 
the function in Eq. f|T2|) also describes the propagation of a charged scalar in an electric 
field, and our previous study demonstrated that such fits can be carried out 13 . 



III. LATTICE DETAILS 

To demonstrate our method for extracting nucleon magnetic moments and electric 
polarizabilities from lattice two-point functions, we have employed an ensemble of anisotropic 
gauge configurations with (2-|- l)-fiavors of dynamical clover fermions 2^, 23|. Our ensemble 
consists of 200 lattices of size x /3 = 20'^ x 128. After an initial 1000 thermalization 
trajectories, the lattices were chosen from an ensemble of 7000 spaced either by 20 or 40 
trajectories to minimize autocorrelations. The lattice spacing in the spatial directions is 



as = 0.123(3) fm [2^, |23(], with a non-perturbatively tuned anisotropy parameter of ^ = 
as /at = 3.5, where at is the temporal lattice spacing. The finer temporal spacing is a crucial 
feature for this study, as it allows us to fit reliably more complicated functional forms for 
the time- dependence of correlation functions. For this ensemble, the renormalized strange 
quark mass is near the physical value, while the renormalized light quark mass leads to a 
pion mass of ^ 390 MeV. 

On each configuration, we compute at least 10 propagators for each of the up, down, 
and strange quarks with random spatial source locations. Multiple inversions were made 
efficient using the EigCG inverter implemented in the Chroma lattice field theory library j2^ . 
Interpolating fields at the source are generated from gauge- covariant ly Gaussian-smeared 
quark fields (25l . 26| on a stout-smeared [23] gauge field in order to optimize the overlap 



onto the ground state in the absence of background fields. Interpolating fields at the sink 
are constructed from local quark fields. Each propagator is located with source time at 
(3^4)src = 0. Randomization of the source time location, while improving the statistical 
sampling, would complicate the determination of charged baryon correlation functions, as 
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their two-point functions are no longer time-translationally invariant. The correlator given 
in Eq. ( ITTj) is generally a function of the sink time-slice and not simply a function of the 
source-sink separation. The full dependence on source time is given in |20|. 

To implement the background field on the lattice, we modify the SU{3) color gauge 
links, Ufj_{x), for each quark flavor by multiplying by the color-singlet Abelian links, ul^\x), 
for the external field, namely 

U,ix) ^ [/^(x)f/f (x)f/f -(x), (13) 

where U^\x) = exp[igA^(x)], with q as the quark electric charge in units of e > 0, and the 
vector potential is chosen to be: A^{x) = (0, 0, — £^''^**a;4, 0), as in Eq. ([1]). The additional 
transverse links are given by 

t/f = exp[tqS'^''/3xs5^^S,,,[s-i]. (14) 



These additional links together with the quantization condition for a torus [28 

= (15) 

ensure that the flux thr oug h every elementary plaquette is qd^^^^^ for down quarks, and 
quS^'^^^ for up quarks [29, 30]. In the quantization condition, q^ = —1/3 is the electric charge 



of the down quark, and n must take on integer values. The use of non-quantized fields have 
been investigated in jlH, where it was found that boundary gradients in the field strength 
can lead to energy shifts in the bulk of the lattice as large as the sought after shifts due 
to the polarizabilities. Notice that is given in lattice units, and the conversion to the 
physical electric field strength S is given by 

£ = {eata,)-'£'^'K (16) 

As the gauge field multiplication in Eq. (fT3l) is carried out on pre-existing gauge configura- 
tions, the sea quarks remain electrically neutral. This approximation is imposed because of 
computational restrictions. 

Using Eq. (fT3l) . we computed propagators for nine values of the field strength, n, 
corresponding to the integer appearing in the quantization condition, Eq. (fT5|) . We use 
n = 0, which corresponds to a vanishing external field, as well as n = ±1, . . ., ±4. On these 
lattices, the expansion parameter governing the deformation of a hadron's pion cloud is given 
by 0: {e£/mlf = 0.18 n^. From the size of this parameter, we anticipate the need to 
include terms beyond quadratic order in the electric field expansion of hadron energies. In 
our analysis, we include terms up to quartic order. Larger lattices with their smaller allowed 
field strengths will be required for better control over systematic uncertainty relating to the 
electric field expansion of observables. 



IV. LATTICE RESULTS 

Nucleon two-point functions were obtained for each source location on a given con- 
figuration. Results for multiple source locations on each configuration were then source 
averaged, yielding lattice correlation functions we denote by: g±{x4,n)i. Here i labels the 
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N: n=0 



0.30 



0.28 




t/a, 



FIG. 1: Effective mass plot for the nucleoli in vanishing electric field (n = 0). We also show the 
effective mass of a two-state fit to the lattice correlation function, where the error band refiects 
only the uncertainty in the ground state energy. The fiat band shows the value of the extracted 
ground state energy with its uncertainty. 

configuration, and ± refers to the boost projection. This procedure was repeated for each 
value of the external field. We have performed multiple differing procedures to analyze the 
data, of which we detail only one method thoroughly in the text. Consistent results were 
obtained from the other procedures. 

To enforce invariance under parity transformations, for which £ — )■ we took the 
geometric mean of correlators calculated at n and —n on each configuration.^ Specifically 
from the set of g±{xi, n)i, we form 



for n > 0. This reduces the nine field values to five, corresponding to the integers n = 
0, . . . , 4. This ensemble of correlation functions was then used to generate 200 bootstrap 
ensembles. For the ensemble averaged correlation functions, we use the same notation but 
without a configuration label, namely 0±(x4,n). Fits to the bootstrapped ensemble are 
performed as described below. 

Fits to correlation functions in vanishing electric field are often guided by effective mass 
plots. Ordinarily one looks for a plateau in the effective mass, Mes{t), to ascertain when 
the excited state contributions have dropped out of the correlator. For the boost projected 
correlation functions, we define two different effective masses 



The situation is quite simple for the case of vanishing electric field; thus we handle this case 
first. 

When the external field vanishes, the two effective masses are identical, Mgjj(t) = 
M^{t) = Meff(t), and the standard analysis applies. In Fig. [H we show the effective mass 
plot for the nucleon in zero electric field. Statistical noise dominates the correlator beyond 
the window of time depicted. As we are limited in statistics, we perform a two-state fit to 
extract the mass of the ground state. The fit function, Q{t,n = 0), has the form 



^ Consistent results were obtained by performing the analysis on the arithmetic mean of correlators. 




(17) 




(18) 




(19) 
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where the parameters Z{0), and M arise from the ground state, while the primed parameters 
account for excited state contamination. We use a correlated chi-squared analysis to fit the 
time dependence of the bootstrap ensemble of correlation functions. As the amplitude 
parameters Z{0), and Z'{0) enter the fit function linearly, we utilize variable projection 
(see [s^] for references) to reduce the number of fit parameters from four down to two. The 
fit to the zero field nucleon correlation function has also been shown in the figure. The fit 
window has been determined by comparing single and double effective masses, see 32N34| for 
details on the latter. The ground state mass we extract from the two-state fit is consistent 
with the high statistics study js^]- 



A. Neutron 



For non-vanishing electric fields, fit functions for the neutron and proton differ consid- 
erably. For the neutron, the fit function is similar in form to the zero-field case, however, 
there are two distinct fit functions corresponding to the boost projectors V±, namely 



Z{n) [1 ± /i^^"(n)^''^*V^] exp -tE{n)^/l - (/ii^"(r2)f '^"/O^ 
+Z'{n) [1 ± /i"^**(n)fi^"/e] exp \-tE'{n)^l - {fi'^^^^{n)£^'''y^y 



(20) 



The unprimed parameters are those of the ground state, while the primed parameters ac- 
count for excited state contributions. Notice both fit functions are identical to Eq. ( IT9l) for 
vanishing electric field. For a fixed non-zero value of the electric field strength, there are 



six fit parameters, three for the ground state: Z{n), E{n), jj, 



latt, 



n] 



and similarly three for 



the excited state contribution.^ To perform the fits, variable projection is again utilized to 
remove the overall amplitudes, Z{n) and Z'{n). This reduces the number of fit parameters 
from six down to four. To determine the remaining four parameters, we perform simultane- 
ous fits to both boost projected correlators for each value of the external field. ^ In principle, 
such fits should take into account correlations between the boost projected correlators. We 
find, however, that the off-diagonal correlations between boost-projected correlators are an 
order of magnitude smaller than the diagonal ones. Thus we treat the boost-projected 
correlators as uncorrelated, and fit them to the function in Eq. (1201) taking into account cor- 
relations in time. In Fig. |2l we show the effective mass plots for the boost projected neutron 
correlation functions. Along with these plots, we show the effective masses resulting from 
the simultaneous fit to both boost projected correlators. Details of the fits to correlation 



^ In principle, the amplitude Z(n) may be different for the differing boost projections, Q+{t, n) and Q-{t, n). 
Any such difference, however, is purely statistical in origin, and a suitable number of measurements should 
produce a common amplitude for the boost-projected correlators within uncertainties. 

^ A simultaneous fit is not required. Alternatively one can separately fit the two boost-projected correlation 
functions, and combine these results to determine energies and magnetic moments. In pursuing this 
alternate procedure, we find that fits to the plus-projected correlation functions are always better than 
fits to the minus-projected correlation functions. For small field strengths, the difference is insignificant; 
however, for the largest field strength, the fit to the minus-projected correlation function is poor. We do 
not presently know the origin of this effect. Using a simultaneous fit to both boost-projected correlators 
mitigates (but does not remove) the problem with the largest field strength. 
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N: n=l 



N: n=l 




t / a, t / a, 

FIG. 2: Effective mass plots for the boost-projected neutron correlation functions. For each value 
of the electric field strength, the curved band shows the result of the simultaneous fit to both 
boost-projected correlation functions using Eq. (j20p . The band accounts for the uncertainty in the 
extracted ground state energy, E{n). The flat band shows the extracted value of E{n) with the 
uncertainty. 

functions, and the extracted parameters are collected in Table HI For the largest value of 
the field strength, £, the confidence of the fit is frankly poor, and so we are careful about 
using this data point. 

The correlation function fits are carried out on each bootstrap ensemble. In partic- 
ular, we arrive at an ensemble of ground state energies and ground state magnetic cou- 
plings for each magnitude of the electric field or equivalently the corresponding integer n. 
These ensembles we generically denote by where i indexes the bootstrap ensemble, 

i = 1, . . . N , and O represents either the ground state energy or magnetic coupling /i. As 
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TABLE I: Summary of fit results for neutron two-point functions using the time window: 5 < t/at < 
28. All quoted values are averages over the bootstrap ensemble, and are given in dimensionless 
lattice units. The conversion to physical units is detailed in the text. For the fits, /d is the 
minimized value for chi-squared per degree of freedom, and 1 — P is the integrated chi-squared 
from the minimum value to infinity. The first half of the table summarizes the time-correlated 
fits to the energies and magnetic couplings in each field using Eq. (j20p . while the second half 
summarizes the field-correlated fits using Eqs. ()2ip and ()22p . The two differing fits to the latter 
are denoted by I and II, and are described in the text. The second uncertainty on polar izabilities 
and magnetic moments is an estimate of the systematic due to the choice of fit window. 



N n 


atE{n) 


^latt 


xVd 


1 - P 





0.206(22) 




0.55 


0.96 


1 


0.2074(16) 


-51(6) 


0.70 


0.93 


2 


0.2142(16) 


-52(3) 


0.91 


0.65 


3 


0.2240(15) 


-50(2) 


1.1 


0.24 


4 


0.2375(15) 


-47(1) 


1.5 


0.02 



N 


atM al" 


xVd 1 - P 


latt 


xVd 1 - P 


I 


0.206(2) 40(9)(2) 


0.3 0.9 


-52(2)(1) 


0.6 0.7 


II 


0.205(2) 42(19)(2) 


0.3 0.9 


-52(3)(1) 


0.7 0.6 



the ensembles of configurations for different field strengths are generated from the same un- 
derlying lattice configurations, correlations between the energies for different field strengths 
will be significant and we account for these. On the bootstrap ensemble of energies and 
magnetic couplings, we perform electric- field correlated fits to the function 0{n), where for 
the case of the ground state energy, 0{n) = E{n), with 



E{n)=M + (S^' 



a 



EEE 



and for the case of the ground state magnetic coupling, 0{n) = fi{n), with 

latt I —latt /'c'latt\2 i —latt /'c'latt\4 
-\- fl^ [C J "I" H-EEE ) ■ 



(21) 



(22) 



With the ensemble average quantitites denoted by 0{n) = minimize the 

correlated chi-squared, namely 



X 



J2 \0{n) - 0{n) 



-1 

n,n' 



C(n') - O(n') 



(23) 



with the field-strength correlation matrix, Cn,n'i given by 

TV 



Cn,n' 



1 



N- 



i=l 



-J2 [0{n)-0,{n)\ [0{n') - 0,{n') 



(24) 



Because all fit parameters enter the fit functions 0{n) linearly, the chi-squared minimiza- 
tion can be done analytically. Fits to the energy function are carried out on the bootstrap 
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N: n 




FIG. 3: Electric field strength dependence of extracted neutron parameters. The two different 
field-correlated fits [I (left panels) and II (right panels)] are described in the text. The bands in 
the plots reflect the total uncertainty. 



ensemble, resulting fit parameters are averaged, and the uncertainties from fitting and boot- 
strapping are added in quadrature. The same is done for the magnetic moment function, as 
defined in Eq. ( 12^ . We find that the best fits result from taking /I'l** = 0, and results quoted 
for the neutron use this constraint. Furthermore we perform two different field-correlated 
fits as follows: (I) a fit to all five field strengths using Eqs. ( 12T|) and fl22|) . (II) the same fit 
function but excluding the largest field strength for which the quality of fit to the correlation 
functions is poor. Finally, to estimate the systematics due to the choice of fit window, we 
performed uncorrelated fits to the electric field dependence of meson energies determined 
on adjacent fit windows. We chose the nine fit windows obtained by varying the start and 
end times by one unit in either direction. On each time window, we determined the elec- 
tric polarizability and magnetic moment. The systematic uncertainty on these observables 
due to the fit window is estimated as the standard deviation of the extracted observables 
over the various adjacent windows. Details of the correlated electric field fits and extracted 
parameters are tabulated in Table [H 

From the extracted parameters, we can investigate the electric field dependence of the 
energies and magnetic moment couplings. This is done in Fig. [3l where we plot the field 
strength dependence of these quantities. The plots, moreover, show the results of the two 
fits (I and II) to the electric field dependence. The values of the extracted parameters are 
consistent with naive expectations and it is useful to convert to physical units. Comparing 
the fit function in Eq. (120]) to the correlator in physical units, Eq. we have 

/i = 2^(»^-^)' ^iatt ^ 0.0313(7) X /i"^", (25) 
with the physical mangetic moment, /i, given in units of nuclear magnetons, /iat — ^ 



2Mn ' 

where Mj^ is the physical mass of the nucleon, and the uncertainty arises from scale setting.^ 



Without a factor of atM/atMN, the magnetic moment would be given in units of lattice nuclear mag- 
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For the magnetic moment of the neutron in units of nuclear magnetons, we thus find 



390 MeV) = -1.63(10)(4)(5) [fi 



N 



We have appended a superscript to reflect that our computation includes only connected 
contributions. The three uncertainties quoted are from: (i) statistics and fitting, (ii) the 
systematic due to the fit window, and (iii) conversion to units of physical nuclear magnetons. 
For (i), we take the largerst value of the uncertainty from the two fits to the field-strength 
dependence (which both gave the same value for /i'*^"). There are additional sources of 
systematic uncertainty that we have not unaccounted for, namely the effects of finite lattice 
spacing and finite lattice volume. 

To convert the lattice electric polarizability to physical units, we compare the fit func- 
tion in lattice units, Eq. fl2T]) . to the energy in physical units, Eq. ([3]), and find 



aE = —atalaf' = 0.0776(58) x af' x 10-^fm^ 

ZTT 



(26) 



where the uncertainty arises from scale setting, and is specifically three times the uncertainty 
in the lattice spacing. For the neutron electric polarizability, we thus find 



390 MeV) = 3.3(1. 5)(2)(3) x lO^^fm^ 



taking the central value and uncertainties from fit II. The notation and sources of uncertainty 
on the electric polarizability are as for the magnetic moment. 

B. Proton 

For the proton, we perform a similar analysis. Fits to proton correlation functions 
are carried out using a two-state fit function. This allows us to remove excited state con- 
tamination. We perform simultaneous time-correlated fits to both boost projected proton 
correlation functions using the fit function 



Zin) 



l±/i^""(n 



latt 



+Z\n) 



1 ± fi' 



(n) 



D t,E(n) 



£^latt 



1 



D t,E'(n? 



latt 



£ 



£ 



latt 



latt 


£^latt 




' i 




£^latt 







£ 



latt 



(27) 



with D{x,E'^,£) as the relativistic propagator function given in Eq. ( IT2|) . As the overall 
amplitudes Z{n) and Z'{n) enter the fit function linearly, we utilize variable projection 
to eliminate them from the simultaneous fits. For each value of the electric field £ (or 
equivalently the integer n), there are then four parameters in the fit: the ground state rest 
energy, E{n), the ground state anomalous magnetic coupling, /i(n), as well the rest energy 



netons, /i'^** = with M as the lattice value of the nucleon mass. With these units, there is no 
uncertainty from scale setting, however, they introduce additional pion mass dependence of the extracted 
moment. 
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FIG. 4: Effective mass plots for the boost-projected proton correlation functions. For each value 
of the electric field strength, the curved upper band shows the result of the simultaneous fit to 
both boost-projected correlation functions using Eq. (j27p . The band accounts for the uncertainty 
in the extracted ground state rest energy, E{£). The curved lower band shows the contribution to 
the fit from only the ground state, while the flat band shows just the extracted value of E{£) with 
its uncertainty. 



and anomalous magnetic coupling for the excited state. Note that we force the proton 
charge to have the value Q = 1. Because we have an improved current, we expect only 
O{o?) differences from the continuum value. In Fig. HJ we show the effective mass plots 
for the boost projected proton correlation functions. Along with these plots, we show the 
effective masses resulting from the simultaneous fit to both boost projected correlators using 
Eq. fl27|) . Details of the fits to proton correlation functions, and the extracted parameters 
from the fits are collected in Table [Tll 
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TABLE II: Summary of fit results for proton two-point functions using the time window: 6 < 
t/at < 28. Tabulated entries are as in Table HJ with the exception that we have denoted the 
magnetic couplings as anomalous using tildes. 



P n 


atE{n) 


^latt 


xVd 


l-P 





0.2052(24) 




0.61 


0.93 


1 


0.2072(17) 


53(6) 


0.65 


0.97 


2 


0.2118(22) 


48(3) 


0.81 


0.80 


3 


0.2198(26) 


46(2) 


0.86 


0.74 


4 


0.2293(29) 


41(1) 


1.5 


0.02 



P atM a"!" x^/d l-P 


rlatt 


xVd 1 - P 


I 0.205(2) 32(13)(1) 0.14 0.98 

II 0.205(2) 31(25)(4) 0.16 0.96 


52(3)(1) 
52(4)(1) 


1.3 0.3 
1.7 0.2 




We perform fits to proton correlations functions on the entire bootstrap ensemble. This 
enables us to form an ensemble of extracted parameters for each value of the field strength, 
In particular, we consider the ensemble of extracted ground state rest energies, {Ei{n)}, 
and ground state anomalous magnetic moments, {fii{n)}. Collectively we denote these 
ensembles by Oi{n), with the ensemble average denoted by 0{n). Electric field correlated 
fits are performed using the fit functions in Eqs. f l2T|) and fl22|) . For the latter it is the electric 
field dependence of the anomalous couplings that is being fit. Furthermore, we extract the 
anomalous magnetic moment using the constraint JI^eee = 0; which results in better fits. 
For both observables, we perform the fit using all the data (fit I), and excluding results for 
the largest field strength (fit II). Results of the fits are collected in Table HTl The electric 
field strength dependence of the extracted rest energies and anomalous magnetic couplings 
are shown in Fig. [51 Also depicted are the field-correlated fits to these quantities. The values 
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of the extracted parameters are again consistent with naive expectations. For the magnetic 
moment of the proton, converting to units of nuclear magnetons using Eq. (125|) yields 



/i™^"(m^ = 390 MeV) = 2.63(13) (1) (4) [^^]. 

We have appended a superscript to reflect that our computation includes only connected 
contributions; we have also added in the Dirac contribution, assuming Q = 1. The three 
uncertainties quoted are as before: (i) statistics and fitting, (ii) the systematic due to the 
fit window, and (iii) conversion to physical units. For (i), we take the central value and 
uncertainty from fit II. For the proton electric polarizability, converting the results of fit II 
to physical units using Eq. yields 



a^"°"^(m^ = 390 MeV) = 2.4(1.9) (3) (2) x lO-^fm^, 

where the last uncertainty arises from scale setting. 

Finally let us compare results for the neutron and proton. Within the uncertainty, 
the connected part of the isoscalar anomalous magnetic moment is consistent with zero. To 
compare with experiment, we require additional contributions from disconnected diagrams 
that we have not determined. The isovector combination of moments, however, does not 
have disconnected contributions due to strong isospin symmetry. For the nucleon isovector 
magnetic moment, we find 

ixy{m^ = 390 MeV) = 4.3(2)(1)(1) [/i^]. 

While this value is smaller than the physical moment, chiral corrections drive the magnetic 
moment downward at masses above the physical value 35, 3^. Studies at additional pion 



masses are necessary to extrapolate to the physical point. The value we obtain, moreover, 
is comparable to values extracted from the current insertion method at similar values of the 
pion mass, see (37| . 

For the electric polarizabilities, our results show both isovector and isoscalar compo- 
nents, however, the latter is the dominant one. This is also seen experimentally and from 
chiral perturbation theory. The smaller isovector component, 

al{m^ = 390 MeV) = -0.9(2.5)(3)(4) x lO-^fm^, 

receives smaller chiral corrections, and is less sensitive to the electric charges of the sea (but 
not independent). While values for the electric polarizabilities of the neutron and proton 
are smaller than experiment, chiral perturbation theory suggests large corrections as one 
nears the chiral limit [38l-41|. Additionally including contributions from sea quark electric 



charges will drive both polarizabilities upwards, as can be seen from partially quenched 
chiral perturbation theory [13]. It will be interesting to carry out simulations at additional 
quark masses and with electrically charged sea quarks to observe this behavior. 



V. CONCLUSION 

Above, we investigate the relativistic propagation of spin-half particles in classical elec- 
tric fields. The presence of magnetic moments affects the behavior of two-point correlation 
functions, and we use this observation to devise a method to determine magnetic moments 
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and electric polarizabilities from lattice QCD simulated in background electric fields. Using 
anisotropic gauge configurations with dynamical clover fermions, we perform such computa- 
tions. In Appendix |Bl we obtain results for the neutron using upolarized lattice correlation 
functions. Such results, however, do not allow one to determine the electric polarizability — 
only a combination of the electric polarizability and the square of the magnetic moment. 
The separation of these terms (which is analogous to accounting for Born-level contributions 
in the Compton scattering amplitude, see Appendix requires treatment of baryon spin, 
which is afforded by studying boost-projected correlators. Our analysis of boost-projected 
lattice correlation functions demonstrates that nucleon magnetic moments and electric po- 
larizabilities can be extracted from lattice calculations in background electric fields. This 
applies to both the neutron and proton. 

There are a number of possible refinements of our computation that would reduce the 
systematic uncertainties. Currently, our calculations are limited to electrically neutral sea 
quarks, and there is a need to remedy this situation. Furthermore, studies on larger volumes 
will not only reduce finite volume effects, but allow the implementation of smaller values of 
the quantized external field strength. Calculations at various values of the quark mass will 
allow for chiral extrapolations to make contact with the physical QCD point. We intend to 
carry out this work in the future. Finally, our approach can be used to study the magnetic 
moments and electric polarizabilities of the remaining members of the baryon octet. 
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Appendix A: Motion-Induced Electric Dipole Moments 

Motion-induced electric dipole moments underly the Born-like terms encountered 
above. These Born couplings are magnetic in origin, and their subtraction is required to 
arrive at physical electric polarizabilities. Born subtractions are carried out when analyzing 
experimental data, calculating polarizabilities from chiral perturbation theory; and, unique 
to this work, determining electric polarizabilities from lattice QCD. Here we remind the 
reader of the physics underlying the Born subtraction. For simplicity, we consider the neu- 
tron in an external electric field. For the proton, there are additional Born couplings to 
the total charge. In momentum space, such terms can be subtracted; whereas in coordinate 
space, these terms must be treated to all orders in the field strength. This is accomplished 
by the relativistic proton propagator in Eq. fITT]) . 

One way to see that the Born term must be subtracted is to consider the electric 
dipole operator, pe- Consider the Minkowski space Hamiltonian derived from the Euclidean 
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effective action for the neutron that appears in Eq. (jl]). The Minkowski space electric field 
we denote by Em- Without subtracting the magnetic moment contribution, the electric 
dipole operator picks up an additional contribution 



Pe 



OH 



dE 



-inaE E 



M 



M 



2M ' 



(Al) 



that does not vanish when the electric field is turned off. This extra contribution is a 
motion-induced effect as can be seen from neutron matrix elements. For a neutron moving 
non-relativistically, 



{pE)^{Niv)\pE\Niv)) 



-AnaE Em + xv + 



(A2) 



we see that the additional term corresponds to a motion-induced dipole moment j43[. In 
the external field, the electric dipole moment contributes to the total energy in the form 



{E) = E 



M 



{pE')dE', 



M 



--AnaEE 



M + 



2M 



V X E 



M 



(A3) 



The second term is readily identified as the interaction energy of the magnetic moment with 
magnetic field seen in the neutron's rest frame: (m) ■ B, with m = and B = v x Em- 
This explains why one sees a motion-induced electric dipole moment in the frame in which 
neutron moves with velocity v. 

Without neutron motion, contributions from the mot ion- induced electric dipole mo- 
ment naively vanish. At second order, however, the mot ion- induced dipole can interact with 
itself via neutron propagation. We must employ a limiting procedure to handle the nucleon 
pole. In a quantum mechanical notation, the shift due to an intermediate-state neutron with 
energy ko = M + ^Mv'^ has the form 



AE 



N{0) 



2M 



K-E 



M 



N{v] 



M -h 



N{v) 



fx 



2M 



K-E 



M 



iV(0) 



(A4) 



The off-diagonal matrix elements in the non-relativistic limit evaluate to half the value of 
the diagonal matrix elements in that limit. We can thus write the energy shift as 



AE = -rh - [v X Em 



- |Mt72 



m - iv X E 



M 



(A5) 



In the limit of zero velocity, a non-vanishing contribution from the motion-induced electric 
dipole moment emerges. This contribution is the same as that derived in Eq. (E]). 



Appendix B: Analysis of Unpolarized Neutron Correlation Functions 

Here we present the analysis of unpolarized neutron correlation functions. On each 
configuration, for each value of the electric field strength, we form the unpolarized, source- 
averaged lattice correlation function, g{x4,n)i. Parity invariance is enforced by taking the 
geometric mean of correlators obtained for a given field value and its negative. Specifically 
we form 

Q{x4,n)i = \/g{x4,n)ig{xi, -n)i, (Bl) 
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N: n=l 



N: n=2 




FIG. 6: Effective mass plots for unpolarized neutron correlation functions. For each value of 
the electric field strength, the curved band shows the result of the fit to unpolarized correlation 
functions using Eq. ()B3p . The band accounts for the uncertainty in the extracted ground-state 
effective energy, Ees(n). The fiat band shows the extracted value of Eeg{n) with the uncertainty. 



for n > 0. This ensemble of unpolarized correlators was then used to generate 200 bootstrap 
ensembles for n = 0, . . . , 4. The average unpolarized correlator is similarly denoted but 
without the subscript referring to configuration number, namely by g(x4,n). The standard 
effective mass is then formed 

and is used to guide spectroscopic analysis of the unpolarized correlators. 

For a given value of the electric field, £, or equivalently the integer n, we extract the 
effective energy, Ees{n) given in Eq. IQ, using a two-state fit function of the form 

g{t, n) = Z{n) exp [-tE,s{n)] + Z\n) exp \-tE'^^[n)\ . (B3) 

Notice without the boost projection, we cannot disentangle the magnetic moment contribu- 
tion. While there are four parameters to fit in Eq. (lB3p . we utilize variable projection to 
eliminate the amplitudes Zin) and Z'{n), leaving just two parameters: the effective energy 
of the ground and excited states. Time-correlated fits are performed, with results shown in 
Fig. O Fit details and extracted parameters are collected in Table IIIII 

Fits are carried out on the entire bootstrap ensemble enabling us to form an ensemble 
of extracted effective energies of the ground state, {Ees,i{n)}. The average of this ensemble 
we denote Ees{n). Using the fit function 

E,s{n) = M + {S^^'^'f + (^''^")^ (B4) 

electric field correlated fits are performed using Eq. fl23l) . We perform two fits: a fit using all 
the data (fit I), and a fit that excludes results for the largest field strength (fit II). Results of 
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TABLE III: Summary of fit results for unpolarized neutron two-point functions using the time 
window: 5 < t/at < 28. All quoted values are averages over the bootstrap ensemble, and are 
given in dimensionless lattice units. The quantity ^^^^ is the pseudo-polarizability defined in 
Eq. ()B5p . For the fits, /d is the minimized chi-squared per degree of freedom, and 1 — P is the 
chi-squared integrated from the minimum value to infinity. The first half of the table summarizes 
the time-correlated fits to the effective energies in each field using Eq. (lB3p . while the second half 
summarizes the field-correlated fits using Eq. (lB4p . The two differing fits to the latter are denoted 
by I and II, and are described in the text. The second uncertainty on the pseudo-polarizability is 
an estimate of the systematic due to the choice of fit window as explained in the text. 



N n 


atEcs{n) 


xVd 


1 - P 





0.2041(23) 


0.50 


0.97 


1 


0.2058(21) 


0.67 


0.88 


2 


0.2082(23) 


0.97 


0.50 


3 


0.2130(21) 


0.70 


0.85 


4 


0.2204(21) 


0.70 


0.85 



N atM x^/d 1 - P 

I 0.205(2) 17(11)(1) 0.53 0.75 

II 0.204(2) 20(25)(3) 0.52 0.72 




FIG. 7: Electric field strength dependence of the neutron effective energy. The two different 
field-correlated fits (I and II) are described in the text, and the bands show the total uncertainty. 

the fits are collected in Table IIIII The coefficient of the term quadratic in the field strength 
is not physically the electric polarizability because it includes Born-level contributions from 
the magnetic moment. For this reason, we call this coefficient the pseudo-polarizability. In 
our choice of lattice units, the pseudo-polarizability is given by 

yilatt ^latt ^t^(,\a.tt\2 (T3^\ 
= «E - ) ' 

with ^ as the anisotropy factor. The systematic due to the fit window is estimated by 
performing uncorrelated fits on the adjacent fit windows obtained by varying the start and 
end times by one unit. 

The extracted value for the pseudo-polarizability is roughly half the size of the electric 



20 



polarizability, see Table HI Because of the sign of the magnetic moment contribution, we 
expect the pseudo-polarizabihty to be less than the electric polarizability. The values of the 
magnetic moment and electric polarizability extracted from boost-projected correlators can 
be used to find a value for the pseudo-polarizability. Using the results of Table [T] in Eq. (IBSp . 
we find 

I : = 17(9)(2), and II : = 19(19)(2), (B6) 

for the two field-correlated fits. These values are concordant with those found in Table IIIII 
from analyzing the unpolarized neutron correlators. 

Finally, we note that unpolarized proton correlation functions in principle allow one 
access to both the magnetic moment and electric polarizability. This can be seen from 
the explicit form for the proton two-point function derived in the effective hadronic the- 
ory. The functional form, however, leads to fits that are challenging to perform. Valuable 
simplifications are afforded by boost projected correlators. 
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